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ABSTRACT 


This thesis addresses the problem of how to locate and staff recruiting stations with 
Active and Reserve recruiters in order to maximize the annual number of recruits. The 
problem is formulated as a nonlinear integer programming problem. The objective 
function for the problem, also referred to as the production function, describes the number 
of recruits obtainable from each zip code and can be estimated via Poisson regression. 
The resulting nonlinear integer programming problem is heuristically solved by 
decomposing decision variables into two sets: one to locate stations and the other to staff 
them with recruiters. Comparisons are made between problems with production functions 


derived from all zip codes and those derived from only zip codes belonging to efficient 


stations as defined in Data Envelopment Analysis. 
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EXECUTIVE SUMMARY 


To support the ongoing drawdown by the Department of Defense, the US Army 
Recruiting Command (USAREC) is in the process of realigning its organizational 
structure for recruiting young men and women to join the Army. Of great concern is the 
question of which stations are to be closed and how to staff the remaining stations with 
recruiters for both the Active and Reserve components. To aid in this decision making, 
this thesis develops an optimization model that takes as its inputs the number of stations 
and the numbers of Active and Reserve recruiters available to a recruiting battalion. Its 
output consists of a list of stations to remain open and the corresponding number of 
Active and Reserve recruiters to staff each of them. 

An integral part of the optimization model is the production function which 
describes the expected number of recruits obtainable from each zip code. This production 
function is not known with certainty and has to be estimated using a statistical technique 
called Poisson regression. To observe the difference in the annual production of recruits 
under the assumption that all recruiters operate in an efficient manner, two types of 
production functions, average and efficient, are considered. The average production 
function is based on data from all zip codes and the efficient one is based on data from 
zip codes belonging to efficient stations. The thesis uses Data Envelopment Analysis to 


determine which stations are efficient. 





To illustrate its utilities, the model was used to locate and staff stations in the 
Albany Battalion with recruiters. It was also observed that a significant number of 
recruits can be obtained if all recruiters are efficient. Although it is optimistic to make 
such an assumption, results from the model with efficient production functions can serve 
as a goal all recruiters should strive to achieve, especially in the current budget 


environment. 
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I. INTRODUCTION 


After forty years of Cold War, when the missions and challenges facing the US 
Armed Forces were clearly defined and easily understood we find ourselves in a period 
of unprecedented change. An increased demand for social and domestic improvement has 
replaced the dissipating threat of the Warsaw Pact. This change in focus brought about 
a corresponding shift of resources, with the Department of Defense being a major target 
for reductions. These reductions affect the number of personnel, the operational funds, 
and the development and acquisition of weapon systems. While the recent number of 
regional conflicts and humanitarian missions indicate that the world remains volatile, the 
reductions will continue. 

The US Army is the most people intensive of all of the Armed Services and 
therefore implementing the personnel drawdown is a point of great concern. To prevent 
the development of a hollow force the drawdown has not been accomplished solely 
through reduced accessions, but rather by making reductions at every level, using a 
variety of incentive and control programs. The budget cuts have been felt through the 
entire force, compelling every unit and organization to become more efficient: being able 


to do more with less. 








A. BACKGROUND 

The drawdown affects the US Army Recruiting Command (USAREC) in several 
ways. USAREC’s primary mission is to recruit young men and women, mainly between 
the ages of 17 and 21, to join the Army. The current downsizing has reduced the 
requirement for new Army recruits from 127,100 in FY92 to 75,000 in FY93. For the 
current fiscal year, as well as next year, USAREC is required to produce 70,000 
enlistments for the Active Army and 46,000 for the Reserves. This reduction has been 
accompanied by smaller recruiting and advertising budgets as well as a smaller recruiting 
force, marked by the elimination of 1,100 recruiters in 1993 alone [Ref. 1]. Meanwhile, 
colleges and other civilian job training institutions have increased their recruiting efforts 
as the population of 17-21 year old individuals is projected to decline by six percent from 
1990 to 1995[Ref. 2]. In addition, today’s emerging weapon technologies demand high 
quality and more capable recruits. These two factors combine to shrink the pool of 
possible recruits for USAREC. Compounding this unfavorable situation is the downward 
shift in the attitude of youths toward a career in the military. During the past three 
years, there has been a 31% decrease in the propensity of young men and women to join 
the military [Ref. 1]. This decline can be attributed to the publicity surrounding the 
continued drawdown, the recent Gulf War and US military involvement in Somalia, and 
other social and economic factors. In order to maintain its competitive advantage over 
other services and civilian organizations in recruiting young men and women, USAREC 


must become as efficient as possible in every facet of its operations. 


B. PROBLEM DEFINITION 

In recruiting, one of the most important resources are recruiters, for they generate 
enlistment contracts for the Army. Therefore, it is important that USAREC provides 
sufficient support for recruiters to perform their duty in the most effective and efficient 
manner possible. In particular, USAREC views recruiting stations as an important 
resource for its recruiters and success in recruiting depends in part on the placement and 
staffing of these stations. A recruiting station provides space for conducting business and 
a homebase for recruiters. Moreover, the presence of a recruiting station also serves as 
an important patriotic reminder in the surrounding community and in some cases attracts 
youths to join the Army. Therefore, USAREC is interested in determining optimal 


locations and staffing levels for its stations. 


C. APPROACH 

This thesis addresses the problem of determining the locations and staffing levels 
of Army recruiting stations in a manner similar to Schwartz [Ref. 3]. The thesis 
formulates the problem as an optimization model with the objective of maximizing the 
total number of yearly enlistments which is statistically estimated from historical data. 
However, this thesis differs from Schwartz in three critical respects. First, Schwartz 
addressed the problem for the Navy Recruiting Command which only recruits for the 
Active component of the Navy. However, USAREC recruits for both the Active and 
Reserve components of the Army and the model in this thesis addresses both of them. 


The Reserve component presents additional complexity in that recruits joining the Army 








Reserve must reside wit : 50 mile radius of his/her assigned Reserve Center. 
addition, recruiters for the Active and Reserve do not necessarily share the same 
recruiting territories. In fact, Reserve recruiters generally must cover more area since 
there are fewer of them to cover the continental United States. Second, Data 
Envelopment Analysis (DEA) [Ref. 4] is used to focus the estimation of the annual 


enlistments on efficient use of resources. Finally, this thesis also employs Poisson 


regression instead of least squares regression to predict the number of yearly enlistments. 


D. THESIS OUTLINE 

In order to allow a thorough understanding of the underlying rationale used in the 
selection of certain techniques and specific explanatory variables, a description of 
USAREC’s organization and current operations is included in Chapter II. Chapter III 
describes and formulates the Army Location Allocation optimization problem. The 
objective function for the problem, also referred to as the production function, describes 
the number of recruits expected from a zip code. Since this function is not known with 
certainty, Chapter IV uses Poisson regression to estimate it. Using DEA to determine 
which recruiting stations are efficient, this chapter concludes with an analysis of two 
different production functions: one using all zip codes and the other using only those zip 
codes that belong to an efficient station. With these production functions, the optimization 
problem in Chapter HI is a nonlinear integer program, a difficult class of problems to 
solve. As an alternative, Chapter V develops a decomposition approach to produce near 


optimal solutions. Chapter V also presents the implementation of the decomposition 


nen 





technique and analyzes the resulting realignment for the Albany Recruiting Battalion. 


Finally, Chapter VI summarizes the thesis and suggest possible areas for future research. 











II. RECRUITING AT USAREC 


This chapter consists of two sections that provide basic information about recruiting 
in the United States Army Recruiting Command (USAREC). The first section provides 
historical information, organization, and structure. The second section describes the 


recruiting operations as they pertain to the problem outlined in Chapter I. 


A. ORGANIZATION AND STRUCTURE 

USAREC is the proponent organization for recruiting young men and women into 
the Active and Reserve Components of the Army and, as such, it is responsible for one 
of the most critical missions of any organization in the Army. It is one of the very few 
organizations that executes its wartime mission on a daily basis. In addition to recruiting 
into the enlisted ranks of the Regular Army (RA) and US Army Reserve (USAR) units, 
USAREC is also responsible for recruiting candidates for other programs such as Officer 
Candidate School (OCS), Warrant Officer Flight Training (WOFT), and Army Nurse 
Corps (ANC). 

In December 1963, a committee commissioned to study all aspects of recruiting for 
the Army found that the organizational structure for recruiting had major inconsistencies 
and was ineffective. As a result, the US Army Recruiting Service was established in 
1964. The organization’s mission also underwent a major revision in the early 1970s 


when the draft ended and an all volunteer force was implemented. This transition brought 





about significant changes in the focus of the entire recruiting process. In 1978, at the 
direction of the Vice Chief of Staff of the Army, USAREC also assumed the mission of 
recruiting for the Army Reserve and became the Total Army’s recruiting organization. 
Currently, USAREC is a field operating agency under the Office of the Deputy Chief of 
Staff for Personnel. In 1993, the Headquarters moved from Fort Sheridan, Illinois to its 
current location at Fort Knox, Kentucky. The current organizational structure of 
USAREC is presented in Figure 1 [Ref. 5]. The different elements of the organization 


will be explained in the subsections below. 





Figure 1. USAREC Organization 





1. Headquarters, USAREC 
Although the mission of USAREC is significantly different from any other 
Army organization, the headquarters and staff operate in much the same manner as any 
major unit. USAREC is commanded by a major general with a deputy commander who 
is a brigadier general and oversees the operations of the Recruiting Brigades. The staff 
is coordinated and led by the Deputy Commander/Chief of Staff and it consists of nine 


major directorates. The organization of the Headquarters is shown in Figure 2 [Ref. 6]. 


HEADQUARTERS 
US ARMY RECRUITING COMMAND 





Figure 2. USAREC Headquarters 


The missions of the directorates involve analyzing, resourcing, and executing 


the current annual recruiting mission. The staff is also involved in the long range 





planning of the entire organization. Of special note is the Program Analysis and 
Evaluation Directorate (PA&E); it is responsible, among many other tasks, for conducting 
analysis that will ensure that all recruiters have the market available to accomplish their 
assigned mission. PA&E provided much of the data used in this thesis and are also the 


intended end user of the methodology presented here. 


2. Recruiting Brigades (Rct Bdes) 
There are currently four Recruiting Brigades dispersed across the country. 
Their locations are shown in Table 1. Each of the brigades is commanded by a Colonel. 
Although the brigade staffs are not as large as the Headquarters’, they still conduct a great 


deal of short term planning and analysis in order to accomplish their specific missions. 


TABLE I. RECRUITING BRIGADE LOCATIONS 












The brigade staff includes two very important branches that do not exist 
separately below the Rct Bde level: the Market Analysis Branch and the ANC Recruiting 
Branch. The Market Analysis Branch dispatches teams to conduct market studies 


(recruiter zone analyses or RZAs) that determine the boundaries of a particular recruiting 











station’s territory. Some of the historical data used in this thesis are drawn from these 
studies. 

The primary purpose of the Rct Bdes is to synchronize the plans and actions 
among the Recruiting Battalions under its control. Under the current alignment, a brigade 


is responsible for eight to thirteen battalions. 


3. Recruiting Battalions (Rct Bns) 

There are currently 40 Rct Bns located in the Continental United States 
(CONUS) and they are predominantly commanded by Lieutenant Colonels. The Rct Bn 
staffs are much smaller than those of the Rct Bdes, and are designed to deal with only 
near term planning and execution. The Rct Bns provide the lowest level dedicated 
planning organization within USAREC. Each Rct Bn controls between four and six 
companies. 

4. Recruiting Companies (Rect Cos) 

There are currently 216 Rct Cos commanded by Captains who have all had 
previous command experience. These command and control organizations are critical due 
to the dispersion of the recruiting stations. An average Rct Co covers an area of 
approximately 10,000 square miles. The Rct Cos represent the link between the policies 
and programs of USAREC and the recruiters at the stations. Their focus is on mission 
accomplishment and on recruiter training. Each Rct Co is assigned four to sixteen 


recruiting stations. 
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5. Recruiting Stations (RS) 

There are currently 1,466 recruiting stations located throughout the United 
States and many of its territories. Typically located in high traffic commercial areas 
(shopping malls and office buildings), they are essentially the liaison between the Army 
and the civilian community. The number of recruiters assigned to a station varies between 
one and nine. A recruiter either recruits for the Active (RA) or Reserve (USAR) 
component, but not both. Some stations also have recruiters whose primary mission is 
to recruit Army nurses. Generally, there is at least one RA recruiter and at most three 
USAR recruiters at every recruiting stations. However, some stations have no USAR 
recruiters. This is because the Reserves have different requirements for its recruits and 
recruiters. First, each recruit must live within a 50 mile radius of his/her assigned 
Reserve Center, where reservists train one weekend of each month. This radius restricts 
the area in which USAR recruiters can recruit. In addition, USAR recruiters are 
sometimes required to recruit for a particular Reserve Center when it has vacancies 
needed to be filled immediately. Finally, RA recruiters mainly recruit individuals with 
no prior military service between 17 and 21 years old whereas USAR recruiters focus on 


a wider population of 17 to 29 year old. 


B. RECRUITING OPERATIONS 
Recruiters operate much like a saleperson selling an Army career to American 
youths. To avoid unnecessary competition and duplication of efforts, USAREC views 


the Continental United States (CONUS) as a collection of zip codes. For Regular Army 
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recruiting, each zip code is assigned to one RA recruiter. A collection of zip codes 
belonging to the same RA recruiter is call a recruiter zone. The recruiting territory of a 
station consists of zones of recruiters who are assigned to the same station. The same 
method also applies to the Reserves. However, because of the previously mentioned 
special requirements, reserve recruiter zones are not generally aligned with the territories 
of the recruiting stations. For areas outside CONUS, the division of zones and territories 
depends on local geographical structure and overseas postal divisions. To simplify our 
presentation, this thesis focuses only on CONUS. 

The Regular Army’s target population of individuals between 17 and 21 years old 
with no prior military experience may be further divided into two major categories: GSA 
and Non-GSA. A GSA recruit is a high school graduate or senior with a category A 
classification that refers to those who score in the upper fifty percentile of the Armed 
Forces Qualification Test (AFQT) test. Last year, 95 percent of 77,600 recruits were high 
school graduates without prior experience and 70 percent scored in the upper 50 
percentile on their AFQT. For the Reserve Army, the target market is larger and includes 
individuals between 17 and 29 years old without regard to prior military experience. 
However, recruits with prior military service are valuable to the Reserve Army, for they 
save training costs and are knowledgeable about current tactics, doctrine, and equipment 
modernizations. These factors are important for keeping Reserve units in synchronization 
with units in the Regular Army. In fact, soldiers separated from the Army are highly 
encouraged to join the Reserve and over 50 percent of recruits that joined the Reserve 


Army in FY93 have prior military service. [Ref. 7] 
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Ii. OPTIMAL ARMY LOCATION AND ALLOCATION PROBLEM 


This chapter presents an optimization problem that determines both the locations for 
recruiting stations and the number of Active and Reserve recruiters for each station. In 
the first section, the problem and its assumptions are stated. The second section provides 
a discussion of prior research related to this type of problem. Finally, the formulation of 


the problem is presented in the last section. 


A. Problem Description 

A set of candidate locations for recruiting stations is assumed known. This is a 
reasonable assumption because downsizing is being considered and the set of candidate 
locations is taken to be the existing station locations. Next, it is assumed that there are 
two production functions, for RA and USAR recruiting, respectively. These functions 
describe the expected number of recruits that can be obtained annually from a given 
zipcode based on (i) demographic and economic factors, (ii) distance to its assigned 
station and (iii) amount of time recruiters (measured, e.g., in man-years) spent recruiting 
in the zipcode. (This recruiting time is also referred to as "recruiter share.") Given this 
information, the problem has four sets of decisions. The first set is to determine which 
candidate stations to open. The second is to assign zipcodes to open stations in order to 
establish the territory of each station. The third is to allocate Active and Reserve 


recruiters to the open stations. Finally, the last set is to decide the recruiter share for each 
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zipcode in a station’s territory. In the optimization problem, these four sets of decisions 


are made to maximize the annual number of Active and Reserve recruits. 


B. Related Research 

Extensive research has been conducted recently on realigning the structure of 
military recruiting organizations. In 1992, Celski [Ref. 8] developed a methodology to 
realign the Army Recruiting Battalions and Companies. In realigning the battalions, his 
model also takes into account state boundaries. When realigning companies within a 
battalion, he assumed that CONUS consists of a collection of counties and his model 
determines which counties belong to which company in an optimal manner. Doll [Ref. 
9] and Schwartz [Ref. 3] addressed problems similar the one described above; Doll’s 
work applied to the Marine Corps and Schwartz’s to the Navy. One key difference 
between our model and those of Doll and Schwartz is the fact that theirs take into account 


only the active component of the respective services. 


C. Problem formulation 


Below is the formulation of the Army Location and Allocation (A-LOCAL) 


problem. 

INDICES: 
S = Candidate Recruiting Station 
Z = Zipcode 
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DATA: 


WA = Weight for Active production function 

WR = Weight for Reserve production function 

NA = Number of available Active recruiters 

NR = Number of available Reserve recruiters 

NS = Number of available recruiting stations 

hie(ar) = Active component production function, where d is the distance 


from zipcode z to its assigned station and r is the recruiter 
share devoted to zipcode z 


g7(d,r) = Reserve component production function 
dz = Distance from zipcode z to station s 
VARIABLES: 

Y, = indicates whether station s is open or closed 

AX7s = indicates whether zipcode z is assigned to station s for 
Active recruiting 

RX; 5 = indicates whether zipcode z is assigned to station s for 
Reserve recruiting 

ASH, = recruiter share devoted to zipcode z for 
Active recruiting 

RSH; = recruiter share devoted to zipcode z for 
Reserve recruiting 

AR, = number of Active recruiters assigned to station s 

RR; = number of Reserve recruiters assigned to station s 
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MAXIMIZE wO2 WA f,(d,,AX,,, ASH,) + WR*g,(d,,RXz5, RSH) 


Vz, s 


SUBJECT TO: 
> Y, = NS 
s 
AX,5 < Ys Vz,s 
Axe s 1 Vz 
Ss 
), AR, = NA 
s 
> (ASH, * AX,,) < AR, Vs 
Zz 
RKo gy S°¥ 5 Vz,s 
> RXz5 <1 Vz 
s 
>> RR, = NR 
s 
)> (RSH, * RX,,) < RR, Vs 
Z 
Ye € {0,13 Vs 
AXeeg G10; Ly KRXs 6 0514 
AR, € {0,1,2...,NA} , RR, € {0,1,2..., NR} 


ASH,20 , RSH, 20 


16 


Vz 


(1) 


(4) 


(5) 


(6) 


(7) 


(8) 


(11) 


(12) 


(13) 





The Active and Reserve objective function weights allow for several possibilities. 
When both WA and WR are one, then the objective is to maximize the expected number 
of Active and Reserve recruits. For other values of WA and WR, the objective function 
represents a weighted combination of the two type of recruits. Setting one of them to 
zero reduces the problem to either maximizing the expected number of Active or Reserve 
recruits. Constraint (1) ensures that NS stations are open. The next four sets of 
constraints pertain to the Active component. Constraints (2) allow zip codes to be 
assigned only to stations that are open. Constraints (3) guarantee that each zipcode is 
assigned to at most one station. Constraint (4) allocates NA recruiters to open stations. 
Constraints (5) apportion recruiter share to each zipcode. Contraints (6) to (9) are for the 
Reserve component and they are analogous to constraints (2) to (5). The remaining 
constraints define which variables are binary, integer and nonnegative. 

The problem as stated above can be applied to the entire CONUS. However, such 
a problem would be too large for many computers. Our implementation in Chapter V 
restricts the problem to the territory of a single battalion. Finally, the A-LOCAL problem 
is a large nonlinear integer programming problem and, therefore, quite difficult to solve. 
A few commercially available software packages, e.g., GAMS/DICOPT [Ref. 10], are 
designed for small to medium size problems. However, none are available to handle a 
problem of this size. Thus, our implementation in Chapter V uses a heuristic approach 


to obtain a good solution to the A-LOCAL problem. 
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IV. FORECASTING RECRUITING PRODUCTION 


One key component of the A-LOCAL problem described in the previous chapter 
is the production functions which estimate the number of recruits for RA and USAR for 
each zip code. In the past, many authors [Ref. 2, 3, and 8] have used standard least 
squares regression to estimate these functions. Least squares regression was the method 
of choice due its wide spread use and its intuitive appeal. Sometimes it provides 
reasonable estimates. Bohn and Schmitz reported that they obtained coefficients of 
determination, R?, between .53 and .60 for their production models using least squares 
regression [Ref. 2]. This low R’ can be in part explained by the fact that least squares 
regression asumes that residuals from the forecasted model are normally distributed; this 
may not be the case in recruiting. In fact, if each individual makes the decision to join 
the Army independently, then the number of recruits from a given zip code has a 
binomial distribution which, in certain limiting cases, can be approximated by either a 
Poisson or normal distribution. However, it is shown below that the Poisson 
approximation is more appropriate for Army data. 

Previous studies have also estimated production functions using data from all zip 
codes. Such an approach does not distinguish efficient recruiters from the inefficient 
ones. This results in production functions that apply to average recruiters -- an “average” 
production function. However, when resources are limited, it is more appropriate to 


estimate the number of recruits that can be obtained by an efficient recruiter -- an 
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"efficient" production function. In fact, an ongoing research project at the Naval 
Postgraduate School is trying to identify factors which will aid in the selection of efficient 
recruiters. Furthermore, data from efficient recruiters may also yield more significant 
relationships between dependent and independent variables. For example, one explanatory 
variable is the distance from a zip code to its assigned recruiting station. The hypothesis 
is that fewer recruits can be obtained from zip codes that are further away from the 
station. For recruiters who do not perform their duty efficiently, distance may not be a 
factor affecting their performance. However, for recruiters who habitually visit potential 
recruits, distance or travel time between the station and zip codes may be a significant 
factor. 

The next section of this chapter describes how to determine efficiency in recruiting 
via Data Envelopement Analysis (DEA). The subsequent section estimates efficient 


production functions based on Poisson regression. 


A. Efficient Recruiting 
In Data Envelopment Analysis (DEA), Charnes, Cooper, and Rhodes [Ref. 4] define 
efficiency for a non-profit organizational unit as the ratio of a weighted sum of outputs 


produced by the unit over a weighted sum of inputs used to produce those outputs, i.e., 


weighted sum of outputs 


Efficiency = oo 
. weighted sum of inputs 
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The weights are scaled so that the maximum value for the ratio is one, representing the 
highest efficiency rating. 

Given the available data, recruiting stations are considered as non-profit 
organizational units for Active recruiting. For our purpose, it would be more precise to 
treat individual recruiters as organizational units. However, there is no accurate data at 
that level of detail. The outputs and inputs for Active recruiting are listed below. The 


list for Reserve recruiting is similar. 


Active Recruiting 
Outputs: 
¢ number of GSA recruits produced by the station 


* number of Non-GSA recruits produced by the station 


Inputs: 
¢ number of RA recruiters at the station 
¢ population of 17-21 year old individuals in the station’s territory 
* number of secondary schools in the station’s territory 
¢ inverse of the area, in square miles, of the station’s territory 
¢ inverse of the average distance from the assigned zip codes to the station 
* average unemployment rate of the assigned zip codes 


¢ average relative military pay in assigned zip codes; defined as the ratio of Army 
base pay to the per capita income 
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In the above lists of inputs, the number of recruiters at each station can be 
changed at the request or the discretion of the station commander. However, inputs 
such as population size and unemployment rate are not under the control of the station 
commander. These inputs are called non-discretionary inputs. [Ref. 11,12] To 
determine the efficiency of recruiting station k, the following optimization, or DEA, 
problem must be solved. 


Data Envelopment Analysis (DEA) Problem 


INDICES: 
S = the Recruiting Station 
i = input for the station 
d = discretionary inputs for the station 
nd = non-discretionary inputs for the station 
O = output for the station 
DATA: 
xX, ; = amount of input used at station s 


Y;,o = amount of output o produced by station s 


VARIABLES: 
u, =the weight given to output o 
v; = the weight given to input i 


> Ug Vk,0 ~ > Vind Xk, nd 
() nd 


MAXIMIZE 
> VaXkd 
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SUBJECT TO: ya » Vice 
oO ni 


<1 Vs 
Y VaXs.a 
d 
u, > 0 Vo 
v,> 0 Vi 


In both the objective function and the constraints, the ratio is slightly different from the 
traditional definition of efficiency used in Charnes et al. [Ref. 4]. In the numerator, the 
weighted sum of outputs is adjusted by the weighted sum of non-discretionary inputs. 
The basic idea in calculating efficiency in DEA is to maximize the efficiency ratio of 
station k subject to constraints that normalize the largest efficiency rating to one. The 
difficulty in the above formulation is in satisfying the requirements that the weights must 
be strictly positive. To handle this difficulty, the technique developed by Springer [Ref. 
11] is used. 

To illustrate the use of DEA, consider the stations of the 1st Recruiting Brigade. 
The majority of the data for these stations came from the Army Territory Assignment 
System (ATAS) database maintained by USAREC. Other demographic information such 
as the unemployment percentage, and the per capita income, are from CACI Marketing 
Systems [Ref. 13]. The problem was solved using the General Algebraic Modeling 
System (GAMS) and a nonlinear program solver called MINOS [Ref. 14]. The results 


are presented in Tables II and If]. To obtain the data presented in the first row of 
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Table II, 32 DEA problems must be solved, i.e. one for each station in Battalion 1A, or 


the Albany Battalion. Data for the other rows in both tables are similarly obtained. 


TABLE II. ACTIVE COMPONENT DEA RESULTS 
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TABLE I. RESERVE COMPONENT DEA RESULTS 





In summary, the purpose of the DEA analysis is to preselect recruiting data that 
represents the work of efficient stations/recruiters. This data will be used in the next 
section to develop efficient production functions. For example, the efficient production 
function for Battalion 1A, the Albany Battalion, will be based on data from the zip codes 
belonging to the 12 efficient stations. On the other hand, the average production function 


will be based on zip codes belonging to all 32 stations in the battalion. 


B. ESTIMATING THE PRODUCTION FUNCTION 

As applied to recruting, production functions give the number of recruits that can 
be obtained from a given zip code. To estimate such functions, this thesis assumes that 
individuals independently decide to join the Army. Also the numbers of recruits from 


distinct zip codes are regarded as independent random variables. Within a single zip 
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code, the accession of individuals can be represented as a sequence of independent 
Bernoulli trials. In this case the number of recruits from each zip code has a binomial 
distribution with parameters n and p. Here, n represents the population of 17 to 21 (or 17 
to 29) year old individuals and p is the probability that a person will join the Army. To 
estimate p, one can assume that it is a function of some explanatory varibles and 
maximize the corresponding likelihood function to obtain the necessary coefficients as in 
logistic regression. Such an approach yields a likelihood function which is nonconcave 
and may produce multiple local optimal solutions. Moreover, it is not clear how to model 
P appropriately under the binomial assumption. An alternate approach is to use the fact 
that, in some cases, the poisson distribution provides a good approximation to the 
binomial distribution. 

When the binomial parameters, n and p, satisfy the conditions: n = 100, p < 0.01 
and np < 20, then the Poisson distribution provides a good estimate to the binomial 
distribution [Ref. 15]. The 1993 data from the Ist Brigade show that the number of 17 - 
21 year olds, with no prior service, varies from zero to 10,801 with an average of 502 
per zip code, and the number of 17 -29 varies from zero to 18,516 with an average of 
1,040 per zip codes. Although the largest number of recruits from a single zip code is 
30 and 27 for RA and USAR recruits, respectively, this number of recruits is only greater 
than 20 in ten of over 7,400 zip codes. Thus, data from this brigade seem to satisfy the 
above conditions. Under the Poisson distribution, A = np represents the expected number 
of recruits from a given zip code. To estimate A for each zip code z, the following model 


based on the Cobb-Douglas production function [Ref. 3] is assumed for RA recruits: 








where 


yy pa efeg Ps uP°AsH,Ps ( = ) Bp 


Zz 


Z 


¢ §, =the number of public secondary schools in the zip code. 


¢ U, =the unemployment rate for the zip code. 
¢ ASH, = the recruiter’s share devoted to zip code z. 


¢ D, =the distance from the zip code’s centroid to its assigned station. 


Zz 


The above model is similar to those described in Bohn and Schmitz [Ref. 2]. Instead of 


least squares regression, the exponents, By, Bs, By, Bsy, and Bp are obtained as follows. 


MAXIMIZE 
LLF =} (-4, + n,1n(a,)) +R 
Z 
SUBJECT TO: 0.01 < Bp < 1.0 
0.01 < By, < 1.0 
where 


vy = ePog Psy Pease Psu ( = ) Bp 


Z 


Zz 


The objective function is simply the log-likelihood function of the Poisson 
distribution [Ref. 16], where R is the remainder term that is constant with respect to B’s. 
The constraints ensure that the resulting model for A, is concave with respect to the 


recruiter share, ASH, and the inverse distance to recruiting stations, (1/Dz). 
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Using the data from the Albany Battalion, Poisson regression is used to estimate the 
number of GSA, NPS and PS recruits. For Active recruiting, the focus is on recruiting 
GSA recruits since they constitute approximatly 75% of total Active recruits. The 
explanatory variables for GSA are as given above. In Reserve recruiting, NPS and PS 
recruits constitiute the total USAR recruits and they are separated here because each 
requires different recruiting tactics/strategies. The additional explanatory variables to 
account for NPS and PS are listed below: 

¢ NRC, = the number of Reserve Centers located within 50 miles of a zip code. 


* RCD, = the average distance to those Reserve Centers within 50 miles of the zip 
code’s centroid. 


As before, data for the explanatory variables for the 1,131 zip codes in the Albany 
Battalion are from ATAS and CACI Marketing Systems. For each type of recruit, two 
productions are estimated using Poisson regression. One is the average production 
function, estimated from all 1,131 zip codes and the other, the efficient production 
function, uses data from zip codes belonging to stations with an efficiency rating of one. 
The resulting exponents are in Table IV. 


TABLE IV. COEFFICIENTS FOR THE ALBANY BATTALION 


LEE 
cea | ase | oso | ok = | on | ore | 


pee Poe ie ff fa ef 
| ps | 301 | oss | a4 | 208 | 582 | 152 [739 | 
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To test how well the data fit the model, chi-square test statistics based on the 
Freeman-Tukey deviates [Ref. 17] were computed and are displayed in Table V along 
with the corresponding p-values. The Freeman-Tukey deviates and the Denominator Free 
goodness of fit test were used to account for the small number of recruits from each zip 
code. Note that the p-values generally lie in a reasonable range (0.1 < p < 0.9) indicating 
that the fit is acceptable for both models, i.e., with efficient and with all zip codes. Later, 
these two models will be compared to assess how they affect the decision to locate 
stations and allocate recruiters. Results for the other battalions of the Ist Recruiting 


Brigade are available in Appendix A. 





TABLE V. SUMMARY OF THE ALBANY BATTALION REGRESSION RESULTS 
Dependent 


Average Production Efficient Production 
function function 
Variable 


aS RET 
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V. IMPLEMENTATION OF THE A-LOCAL MODEL 


From the analysis of Chapter IV, the production functions used in A-LOCAL, f, 


and g,, have the following forms: 





AX 6 <7 
f, (dz ,*AX,,, ASH,) = Pp, (ASH,)° ( a! 
Zs 





pte © 
g, (d,,*RX,,, RSH,) = q, (RSH,) (3 ) 
ZS 


where Pz = eo s,Ps uP 
dz 


and es eRBo ss u Ree NRC, FP anc RCD, "Psa 


where all exponents are estimated via DEA and Poisson regression. With these 
production functions, A-LOCAL is a nonlinear, integer programming problem, a very 
difficult class of optimization problems. This chapter employs a heuristic technique to 
obtain a good solution. 

The heuristic technique is based on the observation that A-LOCAL has two basic 
sets of decision variables: one set locates stations and their territories (i.e., Y,, AX,, and 
RX,, ) and the other allocates recruiters (i.e, AR,, ASH,,, RR,, and RSH,,). These two sets 
of variables are linked mainly by the above production functions. When optimal values 


of one set are known, optimal values for the other set can be determined by solving an 
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independent and smaller subproblem. Thus, the first section below determines the optimal 
locations for stations by assuming that the values of variables in the set for allocating 
recruiters are given. Using these optimal locations for stations, the second section then 
optimally allocates recruiters to the stations. It should be noted that Schwartz also used 
a similar decomposition [Ref. 3]. However, he decomposed his problem into four 
subproblems instead of two. In essence, the technique described below is a streamlined 
and generalized version of Schwartz’, for it uses only two subproblems and applies to 
both Active and Reserve recruiting. Finally, the last section in this chapter describes our 


implementation of the heuristic technique using GAMS. 


A. LOCATING STATIONS 

The heuristic approach for locating recruiting stations first assumes that recruiter 
shares, ASH, and RSH,, have been predetermined in some manner. Note that values of 
ASH, and RSH, implicitly determine the number of Active and Reserve recruiters, AR, and 
RR,, at each station since they must equal the sum of recruiter share assigned to zip codes 
in their territories. Therefore, variables ASH,, RSH,, AR, and RR, can be discarded from 
A-LOCAL, for they become constant and are of no consequence to the problem. This 
reduces A-LOCAL to the following: 


Problem Al: 











MAX > p; ( ASH,)° (= 


Z 


sae >)> q, ( RSH;) a ( as) 


SUBJECT TO: 


Constraints (1), (2), (3), (6), (7), (10), and (11) 
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where ASH, and RSH. 7 are constants representing the predetermined recruiter share. To 
further simplify problem A1, note that, since each zip code z can be assigned to at most 


one station, AX, , can equal one for only a single station, s, and, when y is non-zero, 
AX, Y 
(=) =p (arid Axes for each z. 


Applying similar analysis to RX,, , the objective function of problem Al can be written 
as 
MAX 


Y. p, (ASHz)° (> (2) ax,,) + 0 a, (RSHZ)® (YH (+) ORK g) 
Zz Ss Ags Z Ss des 


Observe that, written in this form, the objective function for problem A1 is linear, thereby 
making it a linear integer program. Moreover, the problem also has a structure similar 
to the well-known uncapacitated plant location problem [Ref. 18]. The distinguishing 
feature of problem A! is that it has two commodities, Active and Reserve. Special 
techniques can be developed to take advantage of this structure, however, these are 
beyond the scope of this thesis. Also, a commercially available solver, the X-System 
[Ref. 19], solves problem Al in a reasonable amount of time (see Section C of this 


chapter). 


ot 








B. ALOCATING RECRUITERS 

Solving problem A1 yields optimal locations and territories for recruiting stations, 
i.e., the optimal values for Y,, AX,, and RX,, are known. Setting these variables to their 
optimal values, i.e., Y,*, AX,,* and RX,,*, reduces and decomposes A-LOCAL into two 


zs ? 


subproblems: one for the Active and the other for the Reserve. 


Problem Bl-a (Active): 


MAXIMIZE Ax oe. * 
>So b, (ASH)? (—,4*) 
Z S ZS 
SUBJECT TO: 
> AR, < NA 
S:Y¥*5=1 
>> (ASH,*AX",,) < AR, Vs: Y*,=1 


Z 


AR,€{0,1,2,...NA} 
ASH, 2 0 Wee 


Problem B1-b (Reserve): 


MAXIMIZE Rx*,. ® 
>>> a, (RsH,)* (— 8 
Z S§ d. 


Zs 


SUBJECT TO: 
>> RR, < NR 


siY¥*s=1 


)> (RSH,*RX*,,) < RR; Vs: Y*,=1 
2 


RR,€ {0,1,2,... NR} 
RSH, 2 0 Vz 
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Consider problem Bl-a. Note that it can be decomposed as follows: 


Master Problem: 


MAXIMIZE > Hs (ARg) 
Si Ye"=1 

SUBJECT TO: > AR, < NA 
S: Y,*=1 


AR, € {0,1,2,....NA} Vs: Y*,=1 


Subproblem: 
MAXIMIZE H,(AR,) = > 6, (ASH,)° 


Z: AXz5*=1 


SUBJECT TO: 


where 





When p is between 0 and 1, the Karush-Kuhn-Tucker (KKT) conditions [Ref. 20], yield 


the following solution to the subproblem: 


G5) 





al 1 
ASH*,=AR,(B,1* / SS 6,7] 
Z: AX*7,5=1 


This solution yields the following objective function value for the Subproblem: 


— ee 
H,(AR,) = > 6, [ARs(B,)7°/ YO (Bz) 1] 
Z: AX"25=1 Zi AX" 7521 
ay 
ns ARS? ( > B, 17°) (1-8) 
Z: AX*z,=1 


Substituting H,(AR;) into the Master problem yields the following problem: 


Problem B2: 
MAXIMIZE Yo, (ARs)? 
S: ¥*,=1 
SUBJECT TO: 
>> AR, < NA 
S:Y"stl1 
AR,€{0,1,2,...NA} VWs: Y*,=1 
where: 


aot. 
6; = ( > p;°*) (1-p) 
Zi AX*z521 
This problem is a nonlinear (integer) allocation problem. However, from the 
Statistical analysis of Chapter IV, p is always between zero and one. Therefore, the 


objective function is concave [Ref. 3] and the problem can be solved optimally by the 


Maximal Marginal Return algorithm [Ref. 21] which is stated below for completeness. 
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Maximal Marginal Return Algorithm 
Step _1: Set AR; = 1 for each s such that Y"; = 1, (Every open station must have at 
least one RA recruiter) 


Step 2: Find the station t with the maximum marginal return of an additional recruiter. 


t = argmax {@, ( (AR,+1)°- (AR,)?) } 
S: ¥,°=1 


Step 3: Set AR, = AR, + 1. 


Step 4: If there are no more recruiters to allocate, stop. Otherwise, return to Step 2. 


problem B1-Reserve can be solved in a similar manner. However, in Step 1, RR; is set 


to zero since there is no requirement for every station to have a Reserve recruiter. 


C. IMPLEMENTATION 

In this section, the A-LOCAL problem is applied to the Albany Recruiting Battalion 
in order to illustrate the technique developed in this chapter. The Albany Battalion is 
the largest battalion in the 1st Recruiting Brigade, containing 1,131 zip codes. The results 
discussed below assume that there are 44 existing stations, of which 36 are to remain 
open. 

From Chapter I'V, the production function for the active component that estimates 


the number of GSA recruits is given below. 
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AX. 0.01 
Fy (d,5*AX, 5, ASH,) = ertig pce U,2 3° ASH,° 68 ( a) 
ZS 


For the Reserve, the production functions consists of two components: one is for recruits 


with prior service and the other is for those without any prior service. 


PS, (d,o*RXo5¢ RSH,) =@3-05,0-06 UP NRCS ROD. * 


0.15 
RX, 


RSH°*7* ( ) 


ZS 


NPS, (dz 5*RXz6, RSH,) = @2-45 Sentra Ce NRCS 49 ROD 


RX 0.09 
RSH,°-®° ( 7 25) 
ZS 


However, in order to apply the MMR algorithm, they are combined into one as follows. 


Cp (d,o*RX, 5, RSH, ) = ( e790 g 0-06 Uz NRO A) RCD, oP RSH 84 


Terres? US NRC 8? RED 8) * ( 


g = 0-15+0.09 
2 


where =0.12 


Note that 5 is simply the average of the distance coefficients for PS and NPS. 

Based on a sample of ten problems, this approach for the Reserves, produces 
answers within 10% of optimality. Considering the fact that the objective function is 
obtained through statistical estimation, solutions within 10% of optimality are judged as 
acceptable. Finally, the values for WA and WR are 0.4906 and 0.5094, respectively; they 


are the fractions of the Active and Reserve recruits for FY93. 
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To locate the 36 stations, problem Al was solved using GAMS with the X-System 
as the integer program solver. For Albany, problem A1 contains 77,539 binary variables 
and 79,756 constraints. It took an average of 7.5 CPU minutes to solve the problem on 
an Amdahl 5995 computer at the Naval Postgraduate school. To allow for an easy 
interface between problem Al and problem B1, the Maximal Marginal Return (MMR) 
Algorithm is also implemented in GAMS. It took GAMS another 33 CPU seconds to 
execute and print solution reports for the MMR algorithm. Recall that the MMR 
algorithm requires no solver. The output from this solution process is shown in Table VI 
below. Columns titled "AUTH AR’ and ’PROP AR’ provide the current and ’optimal’ 
allocations of Active recruiters. Columns titled "AUTH USAR’ and PROP USAR’ have 
similar meaning. Note that stations with zero PROP AR, e.g., station 1A1H, are to be 
closed. The remaining columns give predicted number of recruits in each category: 
GSA, NPS and PS. It should be noted that Schwartz reported that his approach obtains 
solution within 10% of optimality. Since our approach is similar, it is expected that 
similar solution quality is obtained. One method for verifying such a claim involves 
solving the A-LOCAL problem as nonlinear programming problem while ignoring the 
integrality restriction. However, the resulting problems would require an excessive 
amount of computing time due to the large number of variables. In fact, Schwartz 


reported a CPU time of five hours for problems approximately half the size of A-LOCAL. 
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TABLE VI. AN OPTIMAL ALIGNMENT OF THE ALBANY BN 
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D. APPLICATIONS AND RESULTS 
This section demonstrates two possible uses of the A-LOCAL problem. The first 
subsection studies the difference between using efficient and average production functions. 


The other subsection shows how results from solving the A-LOCAL problem with varying 
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number of stations and recruiters can be used to determine the appropriate number of 


stations and recruiters. 


1. Comparison of Efficient and Average Production Functions 

Using efficient and average production functions, the A-LOCAL problem for 
the Albany Battalion was solved using 96 Active and 63 Reserve recruiters and stations 
varying from 18 to 44, Figures 3 to 5 graphically display the results. Figure 3 shows that, 
in terms of combined Active and Reserve recruits, USAREC can obtain an additional 
1,426 recruits per year from the Albany Battalion if the all recruiters are assumed to be 
efficient. Figures 4 and 5 display the number of recruits for Active and Reserve, 
separately. From these figures, efficient recruiters would produce 500 and 926 more 
recruits for Active and Reserve, respectively. Although the assumption that all recruiters 
are efficient is unrealistic, the results obtained using the efficient production functions 
provide USAREC planners with goals against which they can measure their recruiting 


productivity. 
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Figure 3. Total recruits using efficient and average production functions 
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Figure 4, GSA recruits using efficient and average production functions 
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Figure 5. USAR recruits using efficient and average production functions 
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2. Determining the Number of Stations and Recruiters 
Running the A-LOCAL problem with various numbers of stations and 
recruiters produces a series of results that can be used to analyze a wide variety of issues 
involving a single battalion. Presented here are the results of A-LOCAL for five different 
numbers of stations and three different combinations of Active and Reserve recruiters. 
To continue with the idea of efficient recruiting, the results in Table VII are based on 
efficient production functions. 


TABLE VII. RESULTS OF A-LOCAL 


Number of RA | Number of USAR | Number of Stations 





4] 





The curves in Figure 6 show that the number of recruits stops increasing after 
36 stations. This indicates that the Albany Battalion should contain no more than 36 
stations. On the other hand, the differences in the three graphs in Figure 6 also indicate 
that more recruiters mean more recruits. This seems intuitive since it is the recruiters that 
generates recruits not the stations. However, it is also expected that the marginal increase 
due to additional recruiters will level off when the number of recruiters is sufficiently 
large. This is because our choice of production function also models the saturation of the 
market when an excessive number of recruiters is present. 

The curves in Figure 6 also provide information for the appropriate number 
of stations and recruiters. For example, the top most curve indicates that 17 stations, 96 
Active and 63 Reserve recruiters would produce approximately 3,700 total recruits in the 
Albany Battalion. By interpolating between the two top most graphs (see Figure 7), the 
same number of recruits can alternately be obtained with 25 stations and 86 = (96+75)/2 


Active and 56 = (63+50)/2 Reserve recruiters. 
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Figure 6. Results for the Albany Battalion 
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Figure 7. Selecting the Number of Stations and Recruiters by Interpolation 
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VI. CONCLUSIONS 


This thesis addresses the problem of how to improve recruiting by better locating 
stations and staffing them with Active and Reserve recruiters. The problem, which is 
called A-LOCAL, is formulated as a nonlinear integer program with the objective of 
maximizing the total number of recruits. Since the number of recruits is not known with 
certainty, it is modelled with the Cobb-Douglas production function and statistically 
estimated using Poisson regression. To study the effect of efficiency in recruiting, Data 
Envelopment Analysis is used to determine stations which are efficient at recruiting. 
Then, two types of production functions, average and efficient, are estimated. The 
efficient production function is based on data from zip codes belonging to efficient 
stations and the average is based on all zip codes from a battalion. 

Under both types of production functions, the resulting A-LOCAL problem is 
difficult to solve optimally. So, a heuristic procedure is developed to obtain a near 
optimal solution instead. This procedure is based on decomposing the problem into two: 
one is a linear integer program and the other is another nonlinear integer program with 
special structure. The linear integer program deals with locating stations and resembles 
an uncapacitated plant location problem, a problem well-known in the operations research 
literature. The other, the nonlinear integer program, allocates recruiters to open stations. 


This problem has a special structure that allows it to be further decomposed into a master 


and subproblems. The subproblems have closed form solutions, thereby permitting the 
master problem to be solved optimally using the Maximal Marginal Return Algorithm. 

The results from A-LOCAL show that, using the efficient production functions, the 
Albany Battalion can obtain an additional 1,426 recruits or approximately 50 percent more 
than those that can be obtained by using the average functions. In addition, it is also 
demonstrated that 36 stations are sufficient for the Albany Battalion. 

This thesis also identifies the following topics for future research. 

1. As mentioned in Chapter II], A-LOCAL can be applied to CONUS instead of a 
single battalion. However, such an approach would produce an optimization problem too 
large for many existing computers. When applied to CONUS, techniques for 
decomposing A-LOCAL into smaller and more manageable subproblems need to be 
developed. 

2. As formulated, A-LOCAL assumes that stations do not have capacity limitations 
in order to allow for possible expansion of existing stations. Although, it is possible to 
add station capacities to A-LOCAL, it is not clear how the resulting problem can be 
solved in practice. Therefore, exact or approximate solution techniques for handling 


station capacities need to be developed. 
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APPENDIX A POISSON REGRESSION 


A. POISSON REGRESSION IMPLEMENTATION 


$TITLE ACTIVE ARMY REGRESSION MODEL CPT Michael J. Teague 
$STITLE USING POISSON REGRESSION 


* (See Appendice B & C) 


$OFFUPPER OFFS YMLIST OFFSYMXREF INLINECOM{ } MAXCOL 150 
$offlisting 


OPTIONS 
LIMCOL = 0 , LIMROW = 0, SOLPRINT = OFF , DECIMALS = 4 
RESLIM = 100, ITERLIM = 10000, OPTCR =0.1,SEED = 78915; 


A attributes for a zipcode / 
DIST dist to station 
AREA area of zip code 
POP population of 17 to 21 years old in zipcode 
SCHOOLS number of secondary schools in zipcode 
RELPAY _ relative military pay in zipcode in 1990 
UNEMP percent unemployment in zipcode in 1990 
GSA gsa contracts in 1993 
NGSA non gsa contracts in 1993 
REG total ra contracts in 1993 
RECSHR _ share for regular army rec 
EFF efficiency from DEA model 
/ 
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I(A) independent var / DIST, RECSHR, SCHOOLS, UNEMP/ 


D(A) dependent var /GSA, REG/ 


¥ 


SETS 
ZC zip codes / 
$INCLUDE %1.ZID 
RS 


ALIAS (Z,ZC); 


TABLE 
INZIP(ZC,A) information about zipcode 
DIST AREA POP SCHOOLS RELPAY UNEMP GSA 
NGSA REG RECSHR_ EFF 
$INCLUDE %1.DAT 


> 


* BLOCK 1 is without DEA and BLOCK 2 is with DEA 
* ...-USE ONLY ONE AT A TIME---- 


SET ZIP(Z), NZIP(Z); 

* BLOCK 1 divide data into two random groups for cross validation 
ZIP(Z) = YES$(UNIFORM(0,1) GE 1/2); 
NZIP(Z) = NOT ZIP(Z); 


* BLOCK 2 fit the efficient zipcodes and check against the same number 
* — ZIP(Z) = YES$( (INZIP(Z,’EFF’) LE 1) AND (UNIFORM(0,1) GE 1/2) ); 
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* — NZIP(Z) = (NOT ZIP(Z)) * YES$(INZIP(Z,’ EFF’) GE 1); 


* cannot have an independent variable with value of 0 for In transform 
INZIP(Z,D$(INZIP(Z,1T) EQ 0) = 0.01; 


* take In of independent variables for cobb douglas transformation 
PARAMETER LINF(Z,]) natural log of data; 
LINF(Z,D$INZIP(Z,D = LOG(INZIP(Z,D); 


PARAMETER DA(Z) dependent variable; 


VARIABLE 
BO constant term or intercept 
B(I) independent variable coefficients 
LLF log likelihood function 


BY 


POSITIVE VARIABLE 
BO 


2 


B.UP(DIST’) = -0.01; 
B.UP(?RECSHR’) = 1.0; 
B.LOCRECSHR’) = 0.01; 
B.LOCUNEMP’) = 0.0; 
B.LO(?SCHOOLS’) = 0.0; 
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EQUATION 
FUN Log likelihood function for POISSON ; 


-~-EXP(B0)*PRODG, INZIP(ZIP,1I)**B())) 
+DA(ZIP)*(BO+SUM(L,B(D)*LINF(ZIP,])) ) 
); 


MODEL POISSON /FUN/; 
PARAMETER EXPECT(*,*); 


* set up scalers for chi square goodness of fit test 

SCALAR TESTSTAT test statistic for goodness of fit 
DFTESTSTAT denominator free test stat 
CHIPROB _ ’prob the chi-sq > TESTSTA’ 
DFCHIPROB.  ’prob the chi-sq > DFTESTSTAT’ 
CHI10 *chi-squared stat at .10 with d.f. > 40’ 
CHIOS *chi-squared stat at .05 with d.f. > 40’ 
CHI01 *chi-squared stat at .01 with d.f. > 40’ 
DF degree of freedom of the model 


FUN.. LLF =E= SUM(ZIP, 
NORMTEST DRR test percent outside 1.96; 


LOOP(D, 


DA(Z) = INZIP(Z,D); 
SOLVE POISSON USING NLP MAXIMIZING LLF; 











EXPECT(ZIP, ACTUAL’) = DA(ZIP); 
EXPECT(ZIP, EST’) = EXP(BO.L)*PROD(L, INZIP(ZIP,I)**B.L(1) ); 
EXPECT(ZIP,’DRR’) = (SQRT(2+4*DA(ZIP)) 
SOQRT(1+4*EXPECT(ZIP,’EST’)))$DA (ZIP) 

+ (1-SQRT(1+4*EXPECT(ZIP,’EST’)))$(DA(ZIP) EQ 0); 


EXPECT(TOTAL’,’ COUNT’) = CARD(ZIP); 

NORMTEST = SUM(ZIP$(EXPECT(ZIP,’DRR’) LT -1.96 OR 
EXPECT(ZIP,’DRR’) GT 1.96), 1)/CARD(ZIP); 

EXPECT(TOTAL’, ACTUAL’) = SUM(ZIP, EXPECT(ZIP,’ ACTUAL’)); 

EXPECT(TOTAL’,’EST’) = SUM(ZIP, EXPECT(ZIP,’EST’)); 


*DISPLAY EXPECT; 
DISPLAY NORMTEST; 


* Do a CHI Square goodness of fit test 
* DENOM FREE TEST 
DFTESTSTAT = SUM(ZIP, 
SQR( SQRT(DA(ZIP)) + SQRT(1+DA(ZIP)) 
- SQRT(1+4*EXPECT(ZIP,’EST’)) ) 
); 


* CHI SQUARE 
TESTSTAT = SUM(ZIP, SQR( EXPECT(ZIP, ACTUAL’) - EXPECT(ZIP,’EST’))/ 
EXPECT(ZIP,’EST’)): 


* Compare test statistic with chi square 
DF = CARD(ZIP) - (CARD()+1); 


* Calculate CHI Square values using approximation from DEVORE 
CHIiG = DF*POWER( | - 2/(9*DF) + 1.28*SQRT(2/(9*DF)),3); 
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CHI0S = DF*POWER( 1 - 2/(9*DF) + 1.64*SQRT(2/(9*DB)),3); 
CHI01 = DF*POWER( 1 - 2/(9*DF) + 2.33*SQRT(2/(9*DF)),3); 

CHIPROB = 1 - ERRORF((2/(9*DF) - 1 + (TESTSTAT/DF)**(1/3))/SQRT(2/(9*DF))); 
DFCHIPROB = 1 - ERRORF((2/(9*DF) - 1 + 
(DFTESTSTAT/DF)**(1/3))/SQRT(2/(9*DF))): 


DISPLAY LLF.L, BO.L, B.L, TESTSTAT,CHIPROB,DFTESTSTAT,DFCHIPROB, 
DF,CHI10,CHI05,CHI01; 


* Do a CHI Square goodness of fit test on cross validation data set 
* DENOM FREE TEST 
DFTESTSTAT = SUM(NZIP, 
SQR( SQRT(DA(NZIP)) + SQRT(1+DA(NZIP)) 
- SQRT(1+4*EXP(BO.L)*PROD(I, INZIP(NZIP,1**B.L()) ) ) 
); 


* CHI SQUARE 
TESTSTAT = SUM(NZIP, SQR( DA(NZIP) - EXP(BO.L)*PROD(I, 
INZIP(NZIP,])**B.L(D)) )/ 

(EXP(BO.L)*PROD(I, INZIP(NZIP,I)**B.L(1))) 

py 


DF = CARD(NZIP) ; 

CHI10 = DF*POWER( 1 - 2/(9*DF) + 1.28*SQRT(2/(9*DF)),3); 

CHIOS = DF*POWER( 1 - 2/(9*DF) + 1.64*SQRT(2/(9*DF)),3); 

CHIO1 = DF*POWER( 1 - 2/(9*DF) + 2.33*SQRT(2/(9*DF)),3); 

CHIPROB = 1 - ERRORF((2/(9*DF) - 1 + (TESTSTAT/DF)**(1/3))/SQRT(2/(9*DF))); 
DFCHIPROB = 1 - ERRORF((2/(9*#DF) - 1 + 
(DFTESTSTAT/DF)**(1/3))/SQRT(2/(9*DF))); 


DISPLAY TESTSTAT,CHIPROB,DFTESTSTAT,DFCHIPROB, DF,CHI10,CHI05,CHIO1; 
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1. Poisson Results of Active Battalions with DEA 


Ce = [emi [oe [om [roe [or 
1.4824 .0624 3039 .6758 188.2 8890 
2.2768 .3409 .0801 .0782 3128 169.3 .1216 
1.7593 0431 2151 .0838 .6752 174.2 .6272 
2.0730 .0346 poo | 0454 .7351 8888 


Pie [ow Prose [0 [ase [ove [as [ae | ao 
To [ose oar [0 saat [oe [is Pass [0 
Pe [ese [ase [0 [me [or [ase a [oe 
Por [ose Pow [ow | 0 [ow] ao | as | mo 
Pos [os [oer [oe [0 [oon [see [ise [ma 
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2. Poisson Results for Active Battalions without DEA 
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3. Poisson Results for Reserve Battalions with DEA 
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4. Poisson Results for Reserve Battalions without DEA 
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APPENDIX B A-LOCAL GAMS PROGRAM 


$TITLE ARMY LOCATION ALLOCATION CPT Michael J. Teague 
$STITLE  Realign Recruiting Stations, Active and USAR Recruiters for a Bn 

* 

*__-_.------ GAMS AND DOLLAR CONTROL OPTIONS ------------------------------- 

* (SEE APPENDICE B & C) 


$OFFUPPER OFFS YMLIST OFFS YMXREF OFFUELLIST INLINECOM{ } MAXCOL 
130 


$SOFFLISTING 
OPTIONS 
LIMCOL=0 , LIMROW=0 _ , SOLPRINT = OFF, DECIMALS = 4 
RESLIM = 90000, ITERLIM=9999999, OPTCR =0.05 , SEED = 3141; 
*----------- DATA AND SETS-----------------------------------------7--9- "> 
SETS 
A attributes / 
ZIPX x-coordinate of zipcode centroid 
ZIPY y-coordinate of zipcode centroid 


POP1721 active target population of 17 to 21 year olds 
POP1729 reserve target population of 17 to 29 year olds 
SCHOOLS number of secondary schools within the zipcode 
UNEMP unemployment percentage of zipcode in 1990 


NRCSO number of reserve centers within 50 miles of zipcode 
RCDIST average distance to reserve centers within 50 miles 
ASHR share of active recruiters assigned to the zipcode 
RSHR share of reserve recruiters assigned to the zipcode 
DIST distance from zipcode to assigned stations 

GSA number of active GSA contracts in 1993 

NPS number of reserve NonPrior Service contracts in 1993 
/ 


L location / RSID, XCOORD, YCOORD, OPRA, OPAGR / 


I(A) independent variables 
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/ SCHOOLS, UNEMP, ASHR, RSHR, DIST, NRC50, RCDIST, POP1721, 


POP1729 / 


? 


SETS 
ZC zipcodes / 
$INCLUDE ZIPBNIA ZID 
/ 
RC _ reserve centers / 
$INCLUDE RESCTR STA 
/ 
S station rsids / 
$INCLUDE RSIDBNIA STA 
oe 


ALIAS (Z,ZC); 


TABLE 
INZIP(ZC,A) information about zipcode 
ZIPX ZIPY POP1721 POP1729 
UNEMP GSA NPS 
$INCLUDE ZIPBN1A DAT 


’ 


TABLE 
INSTA(S,L) information about recruiting stations 
XCOORD YCOORD OPRA OPAGR 
$INCLUDE RSIDBN1A DAT 


> 


TABLE 
LOCRC(RC,L) location of reserve centers 
XCOORD YCOORD 
$INCLUDE RESCTR DAT 


? 


SCALAR BO intercept for active production function 
RBO intercept for reserve production function 
PBO 


a 





SCHOOLS 





NOPRA total number of active recruiters available to BN 
NOPAGR total number of reserve recruiters available to BN 
NS total number of recruiting stations available to BN 
Wil weight for active production from ATAS 93 

W2 weight for reserve production 

RAD maximum distance from zip to RS_ helps solvability; 


PARAMETER 
BQ) _ coefficients of active production function 
PB(]) 
RB(I)_ coefficients of reserve production function _ ; 


*Scalars 
BO = 1.4824 ; 
RBO = 2.4508 ; 


RAD = 150 : 


*Parameters 

BC SCHOOLS’) = 0.0624 ; 
BCUNEMP’) = 0.3039 ; 
BC ASHR’) = 0.6758 — ; 
B(DIST’) = -0.01 : 
RB(’?SCHOOLS’) = 0.1521 ; 
RBCUNEMP’) = 0.2022 ; 
RBCRSHR’) = 0.6006 ; 
RBCDIST’) = -0.0938 ; 
RBCNRCS50’) = 0.4883 ; 
RB(RCDIST’) = -0.3442 ; 
PBC SCHOOLS’) = 0.0581 ; 
PBCUNEMP’) = 0.4136 ; 
PB(RSHR’) = 0.7390. ; 
PB(DIST’) = -0.1517 ; 
PBCVNRCS50’) = 0.2081 ; 
PBC RCDIST?’) = -0.5825 ; 
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PARAMETER 
DR(ZC,RC) distance from zip to reserve center; 


DR(Z,RC) = 69.171*SQRT(POWER(COS(3.14159* 
(INZIP(Z,’ZIPY’)+LOCRC(RC,’ YCOORD’)) /360)* 
(CINZIP(Z,’ZIPX’)-LOCRC(RC,’XCOORD’)) ), 2) 

+ POWER((INZIP(Z,’ZIPY’)-LOCRC(RC,’ YCOORD’)), 2)); 


DR(Z,RC)$(DR(Z,RC) LT .5) = 0.5; 


SET ZS5(ZC) set of zip codes which are within 50 mi of a reserve center; 
Z5(Z) = YES$(SUM(RC$(DR(Z,RC) LE 50), 1) gt 0); 


INZIP(ZC,’ NRC50’) = SUM(RC$(DR(ZC,RC) LE 50), 1); 


INZIP(ZC,’ RCDIST’ )$INZIP(ZC,’NRC50’) = 
SUM(RC$(DR(ZC,RC) LE 50), DR(ZC,RC)) 
/JINZIP(ZC,’NRCSO0’); 


PARAMETER 
D(ZC,S) distance from zip to recruiting station; 


D(Z,S) = 69.171*SQRT(POWER(COS(3.14159* 
(INZIP(Z,’ZIPY’)+INSTA(S, YCOORD’)) /360)* 
((INZIP(Z,’ZIPX’)-INSTA(S, XCOORD’)) ), 2) 

+ POWER((INZIP(Z,’ZIPY’)-INSTA(S,’YCOORD’)), 2)); 


D(Z,S)$(D(Z,S) LT .5) = 0.5; 


* cannot have an independent variable with value of 0 
INZIP(Z,D$(INZIP(Z,D) EQ 0) = 0.01; 


PARAMETER 
C(Z) constant terms for active production function for each zip 
K(Z) constant terms for reserve production function for each zip 
P(Z) 
CHAT(Z) changed C(Z) based on first approximation 
KHAT(Z) changed K(Z) based on first approximation 
PHAT(Z) CHANGED P(Z) BASED ON FIRST APPROXIMATION 
YCOV(Z) ’=1 if zip is within RAD’; 








C(Z) = EXP(BO)*(INZIP(Z,’SCHOOLS’)**B(SCHOOLS’)) 
*(INZIP(Z,, UNEMP’)**BC UNEMP’)) ; 


K(Z) = EXP(RBO)*(INZIP(Z,’ SCHOOLS’)**RB( SCHOOLS’)) 
*(INZIP(Z, UNEMP’)**RB(? UNEMP’)) 
*(INZIP(Z, NRCS0’)**RBC NRCS50’)) 
*(INZIP(Z,’ RCDIST’)**RB( RCDIST’)); 


P(Z) = EXP(PBO)*(INZIP(Z,’ SCHOOLS ’)**PB( SCHOOLS’)) 
*(INZIP(Z,’ UNEMP’)**PB(? UNEMP’)) 
*(INZIP(Z,’ NRC50’)**PBC NRCS50’)) 

*(INZIP(Z, RCDIST’)**PB(RCDIST’)); 


CHAT(Z) = C(Z)*( (NOPRA*INZIP(Z,’POP1721’)/ 
SUM(ZC, INZIP(ZC,’POP1721’)))**BCASHR’)) ; 


KHAT(Z) = K(Z)*( (NOPAGR*INZIP(Z, POP1729’)/ 
SUM(ZC, INZIP(ZC,’POP1729’)))**RBCRSHR’)) ; 


PHAT(Z) = P(Z)*( (NOPAGR#*INZIP(Z,’POP1729’)/ 
SUM(ZC, INZIP(ZC,’POP1729’)))**PBCRSHR’)) ; 


YCOV(Z) = 1$(SUM(S$(D(Z,S) LE RAD), 1) GT 0); 


VARIABLE 
CONTR _ total number of contracts from this BN 


> 


*POSITIVE VARIABLE 


BINARY VARIABLE 
Y(S) open or close station s 
X(Z,S) assign zip to station for active recruiting 
RX(Z,S) assign zip to station for reserve recruiting 


> 


EQUATIONS 
APPROX obj function for linear approx to total contracts 
TOTSTA limit the number of open stations 
AZIP(Z,S) only assign a zip to an open station for active 
RZIP(Z,S) only assign a zip to an open station for reserve 
ACEACH(Z) _ assign a zip to one and only one station for active 
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RESEACH(Z) assign a zip to one and only one station for reserve 


> 


APPROX.. CONTR =E= W1*SUM(Z$YCOV(Z), 
CHAT(Z)*SUM(S$(D(Z,S) LE RAD), 
— (X(Z,S)*(D(Z,S)**BCDIST’)) ) )) + 
W2*SUM(Z5$YCOV(Z5), 
KHAT(Z5)*SUM(S$(D(Z5,S) LE RAD), 
(RX(Z5,S)*(D(Z5,S)**RB(’DIST”)) ) ))+ 
W2*SUM(Z5$YCOV(Z5), 
PHAT(Z5)*SUM(S$(D(Z5,S) LE RAD), 
(RX(Z5,S)*(D(Z5,S)**PB(DIST”)) ) )); 


TOTSTA.. SUM(S, Y(S)) =E= NS; 

AZIP(Z,S)$(D(Z,S) LE RAD).. X(Z,S) =L= Y(S); 
RZIP(Z5,S)$(D(Z5,S) LE RAD).. RX(Z5,S) =L= Y(S); 
ACEACH(Z)$YCOV(Z).. SUM(S$(D(Z,S) LE RAD), X(Z,S)) =L= 1; 


RESEACH(Z5)$YCOV(Z5).. SUM(S$(D(Z5,S) LE RAD), RX(Z5,S)) =L= 1; 


MODEL STATIONS /ALL/; 


$BATINCLUDE ’XSOPTION INC A’ STATIONS GEN 400 200 200 1 * * 
SOLVE STATIONS USING MIP MAXIMIZING CONTR; 
SCALAR OBJ1, OBJ2, OBJ3; 
OBJ1 = SUM(Z$YCOV(Z), 
CHAT(Z)*SUM(S$(D(Z,S) LE RAD), 
(X.L(Z,S)*(D(Z,S)**BCDIST’)) ) )) ; 
OBJ2 = SUM(Z5$YCOV(ZS), 
KHAT(Z5)*SUM(S$(D(Z5,S) LE RAD), 
(RX.L(Z5,S)*(D(Z5,S)**RBCDIST’)) ) )); 
OBJ3 = SUM(Z5$Y COV(Z5S), 
PHAT(Z5)*SUM(S$(D(Z5,S) LE RAD), 
(RX.L(Z5,S)*(D(Z5,S)**PB(DIST’)) ) )); 
DISPLAY OBJ1, OBJ2, OBJ3; 


* Uses the Max Marginal Return heuristic 


SET RCNT loop indice for number of recruiters assigned 
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/1*500/ 


> 


SCALAR 
NRCT counter for number of recruiters assigned 
MXBFIT the max marginal benefit for a single iteration 
MXFOUND counter to id when the sta with MXBFIT is reached 
BETA coefficient for recruiter share 


ES 


PARAMETER 
MRBFIT(S) marginal benefit of assigning an add recruiter to sta 
STAREP(*,*) report for results 


b 


PARAMETER 
DELTA(S) approx active contracts for sta without ASHR 
RA(S) number of active recruiters to assign to station 


4 


BETA = BC ASHR’); 
DELTA(S)$(Y.L(S) EQ 1) = ( SUM(Z$(X.L(Z,S) EQ 1), 
(C(Z)*(D(Z,S)**B( DIST’)))**(1/(1-BETA)) )**(1-BETA) 

)s 

* Initialize values based on at least one recruiter to each open station 

RA(S) = 1$Y.L(S); 

MRBFIT(S) = DELTA(S)*(2**BETA - 1**BETA)$Y.L(S); 

NRCT = SUM(S, RA(S)); 

LOOP ( RCNT$(NRCT LT NOPRA), 


MXBFIT = SMAX(S$Y.L(S), MRBFIT(S)); 
MXFOUND = 0; 


LOOP ( S$( Y.L(S) EQ 1 AND MXFOUND EQ 0), 
IF ( MRBFIT(S) EQ MXBFIT, 
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RA(S) = RA(S) + 1; 
NRCT = NRCT + 1; 
MRBFIT(S) = DELTA(S)*( (RA(S)+1)**BETA - 
. RA(S)**BETA _); 
MXFOUND = 1; 
); {endif} 
); {end loop} 
); {end loop} 


**** REPORT RESULTS **** 

STAREP(S,’OLD OPRA’) = INSTA(S,’OPRA’); 
STAREP(S,’NEW OPRA’) = RA(S); 

STAREP(S,’GSA’) = (DELTA(S)*(RA(S)**BETA))$Y.L(S); 
STAREP(S,’NZIPS’) = SUM(Z, X.L(Z,S) ); 
STAREP(TOTAL’,’OLD OPRA’) = SUM(S, INSTA(S,’OPRA’)); 
STAREP(’TOTAL’,’NEW OPRA’) = SUM(S, RA(S)); 
STAREP(TOTAL’,’GSA’) = SUM(S, STAREP(S,’GSA’)); 


PARAMETER 
RDELTA(S) APPROX RESERVE CONTRACTS FOR STA WITHOUT RSHR 
RDEL1(S) APPROX RESERVE CONTRACTS FOR STA WITHOUT RSHR 
RDEL2(S) APPROX RESERVE CONTRACTS FOR STA WITHOUT RSHR 
USAR(S) NUMBER OF RESERVE RECRUITERS TO ASSIGN TO STATION 


* Uses the Max Marginal Return heuristic 
BETA = RB(’RSHR’); 


RDEL1(S)$(Y.L(S) EQ 1) = SUM(Z5$(RX.L(Z5,S) EQ 1), 
(K(Z5)*(D(Z5,S)**RB(DIST’)))**(1/(1-BETA)))**(1-BETA); 


BETA = PB(RSHR’); 


RDEL2(S)$(Y.L(S) EQ 1) = SUM(Z5$(RX.L(Z5,S) EQ 1), 
(P(Z5)*(D(Z5,S)**PB( DIST’ )))**(1/(1-BETA)))**(1-BETA); 


BETA = (RBCRSHR’) + PBCRSHR’))/2; 
RDELTA(S)$(Y.L(S) EQ 1) = SUM(Z5$(RX.L(Z5,S) EQ 1), 
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(K(Z5)*(D(Z5,S8)**RB(’DIST’))+P(Z5)*(D(Z5,S)**PB(DIST’ )) 
)**(1/(1-BETA)))**(1-BETA) 


> 


DISPLAY RDELTA, RDEL1, RDEL2; 

* Initialize values based on no reqd min recruiters at each open station 
USAR(S) = 0; 
MRBFIT(S) = RDELTA(S); 
NRCT =0; 


LOOP ( RCNT$(NRCT LT NOPAGR), 


MXBFIT = SMAX(S$(Y.L(S) EQ 1 AND USAR(S) LT 3), MRBFIT(S)); 
MXFOUND = 0; 


LOOP ( S$( Y.L(S) EQ 1 AND MXFOUND EQ 0 AND USAR(S) LT 3), 
IF (MRBFIT(S) EQ MXBFIT, 
USAR(S) = USAR(S) + 1; 
NRCT = NRCT + 1; 
MRBFIT(S) = RDELTA(S)*( (USAR(S)+1)**BETA - 
USAR(S)**BETA _ ); 
MXFOUND = 1; 
); {endif} 
); {end loop} 
); {end loop} 


4+ REPORT RESULTS **** 


STAREP(S,’OLD OPAGR’) = INSTA(S,’OPAGR’); 
STAREP(S,’NEW OPAGR’) = USAR(S); 

STAREP(S,’NPS’) = (RDEL1(S)*(USAR(S)**RB(RSHR’)))$Y.L(S); 
STAREP(S,’PS’) = (RDEL2(S)*(USAR(S)**PB(RSHR’)))$Y.L(S); 
STAREP(S,’NRZIPS’) = SUM(Z, RX.L(Z,S) ); 
STAREP(’TOTAL’,’OLD OPAGR’) = SUM(S, INSTA(S,’OPAGR’)); 
STAREP(’TOTAL’,’ NEW OPAGR’) = SUM(S, USAR(S)); 
STAREP(’TOTAL’,’NPS’) = SUM(S, STAREP(S,’NPS’)); 
STAREP(’TOTAL’,’PS’) = SUM(S, STAREP(S,’PS’)); 


FILE SUMMARY /’%4 REPORT’/; 
PUT SUMMARY 
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PUT’ + BATTALION tA. 8° 8 7/5 


PUT’ ts 
PUT’ AUTH PROP PRED AUTH PROP PRED PRED’/; 
PUT’STATION OPRA OPRA GSA OPAGR OPAGR NPS __PS’/; 
LOOP(S, 

PUT S.TL:10, 


STAREP(S,’OLD OPRA’):4:0, STAREP(S,’NEW OPRA’):10:0, 
STAREP(S,’GSA’):10:2, 
STAREP(S,’OLD OPAGR’): 8:0, STAREP(S,’NEW OPAGR’):8:0, 
STAREP(S,’NPS’):10:2, STAREP(S,’PS’):9:2 /; 
* PUT STAREP(S,’NZIPS’):15:0, STAREP(S,’NRZIPS’):10:0/; 
); 
PUT ’TOTAL’:10, 
STAREP(’TOTAL’,’OLD OPRA’):4:0, STAREP(’ TOTAL’, NEW OPRA’):10:0, 
STAREP(TOTAL’,’GSA’):10:2, 
STAREP(’TOTAL’,’OLD OPAGR’): 8:0,STAREP(’TOTAL’,’ NEW OPAGR’):8:0, 
STAREP(’TOTAL’,’NPS’):10:2, 
STAREP(TOTAL’,’PS’):9:2; 
PUT // "STATION CONSTRAINT : " NS; 
PUT / "OPRA CONSTRAINT :" NOPRA; 
PUT / "OPAGR CONSTRAINT : " NOPAGR; 
PUTCLOSE SUMMARY; 
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